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ABSTRACT 

Owing to the frequency and reproducibility of its outbursts, the black-hole candi¬ 
date GX 339-4 has become the standard against which the outbursts of other black- 
hole candidate are matched up. Here we present the first systematic study of the 
evolution of the X-ray lags of the broad-band variability component (0.008 — 5 Hz) 
in GX 339—4 as a function of the position of the source in the hardness-intensity 
diagram. The hard photons always lag the soft ones, consistent with previous results. 
In the low-hard state the lags correlate with X-ray intensity, and as the source starts 
the transition to the intermediate/soft states, the lags first increase faster, and then 
appear to reach a maximum, although the exact evolution depends on the outburst 
and the energy band used to calculate the lags. The time of the maximum of the lags 
appears to coincide with a sudden drop of the Optical/NIR flux, the fractional RMS 
amplitude of the broadband component in the power spectrum, and the appearance 
of a thermal component in the X-ray spectra, strongly suggesting that the lags can 
be very useful to understand the physical changes that GX 339-4 undergoes during 
an outburst. We find strong evidence for a connection between the evolution of the 
cut-off energy of the hard component in the energy spectrum and the phase lags, 
suggesting that the average magnitude of the lags is correlated with the properties 
of the corona/jet rather than those of the disc. Finally, we show that the lags in GX 
339-4 evolve in a similar manner to those of the black-hole candidate Gygnus X-1, 
suggesting similar phenomena could be observable in other black-hole systems. 

Key words: accretion, accretion discs, black hole physics. X-rays: binaries. 


1 INTRODUCTION 


The X-ray spectrum of black-hole candidate (BHC) systems 
in low-mass X-ray binaries (L MXBs) can be decomposed 
into t wo main components: (e.g., iTanakall 19891 : Ivan der Kli^ 
[Tiii): a hard component, usually fitted w ith a power law 
with p hoton index in t he range 1.4-2.5 (e.g., Thome &: Pried 
I 1975 I : ISunvaev et al.l Il99ll : iGrebenev et alT 19931 ). and a 
soft, thermal, blackbo dy-like component with temperature 
kT ~0.5 —2keV (e.g.. lMitsuda et aI.lll98J : lMivamc)to et al.l 
Il993h . The hard component is usually attributed to a corona 
of energetic electrons, and the state where it dominates is 
known as the low-hard state (LHS). The soft component 
is attributed to thermal emission fro m an optically thick 
but g eometrically-thin accretion disc (IShakura fc SunvaevI 
when the soft component dominates the emission, 
the source is in the so-called high-soft state (HSS). 


* E-mail: d.altamirano@soton.ac.uk 


Any description of the spectral evolution during an 
outburst based on spectral fittings is of course model- 
dependent, therefore many authors have opted for using 
X-ray colours (also known as hardness); these colours are 
defined as the ratio between count rates in different bands, 
and have proven to be very useful to understand the spec¬ 
tral evolution of outbursts in terms of the dominating com¬ 
ponent of the spectrum. For BHCs, the hardness-intensity 
diagram (HID) is generally used to characterize the spectral 
properties of a source. In what is now known as “canonical 
outburst”, th e source evolves in a roughly square pa ttern in 
the HID (e.g.. lHoman et ^l200ll : lBelloni et al.ll2005l . among 
many others). First, the intensity of the source increases 
by a large factor while the hardness slightly decreases. The 
source is in the LHS and traces a roughly vertical line in the 
HID. At some point the source “turns the corner”: the in¬ 
tensity stops increasing and the source spectrum softens at 
more or less constant intensity towards the HSS. In between 
the LHS and HSS the source passes through the hard- and 
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soft-intermediate states (see, e.g., iHoman fc Bellonill2005ll . 
showing complex behaviour including sometimes large flares 
in intensity. After the source reaches the HSS, the intensity 
usually decreases at approximately constant hardness, with 
some excursions to harder or softer states; as the outburst 
progresses, the source returns to the LHS at approximately 
constant intensity. The hardening of the spectra continues 
until the hard colour returns to a value similar to those ob¬ 
served at the beginning of the outburst. At this point the 
hardening stops, the intensity drops and the source returns 
to quiescence. 


While the detailed evolution of BHC outbursts can be 
more complex than what we have just described, change 
slightly from source to so urce, or in differen t outbursts of 
the s a me source (see, e.g. , lHomnnet_alJj200iyBellom_e^_aL 
200E; Homan fc BelloiM l2005l : Remillard fc McClintock 


200f; Ivan der KlisI ~06l . and references therein), the loop 
in the HID is usually recognizable. This “ q-shaped” loop is 


q-sn 

usua lly discussed in terms of hysteresis (e.g.lMivanioto_^£_aI 


_1995) ■ and has been known for a long time fe.g. lHoman et al.i 
2001 : iMaccarone fc CoDPill2003ll . We note however, that al¬ 


though we describe this hysteresis effect only in the case of 
BHs in LMXBs, similar loops can be seen in HIDs and/or 
colour-colour diagrams of other accreti ng objects (e.g., neu¬ 
tron s tars and Dwarf Nova e; see e.g., IMaccarone fc Coppil 
I2OO3I : iKording et al.l l2008h . therefore pointing at a simi¬ 
lar physical origin. The accretion physics behind the “q” 
shape in the HID of BHs, and what produces the cor¬ 
ners is n ot fully understood (but see , e.g., proposed sce- 
narios bviMever-Hofmei ster et al.ll2005l:lPetrucci et al.l l200i : 


iBegelman fc Armitagd 2014 : Kvlafis fc BellonT 2015 . and 

references therein); however, it is probably related to sig- 
nihcant changes in the geometrical/structural conhguration 
of the system triggered by changes in mass accretion rate 
(M). For example, there is clear evidence that the ejection of 
relativistic jets takes place during som e of the transitions be¬ 
tween the states described above (e.g. , Vadawale et al.ll2003l : 
iFender fc Belloiiill2004l : ICorbel et al’lf^Odl : lBellonill2007l l~ 


To fully describe the evolution of a black-hole out¬ 
burst it is essential to consider the X-ray (time) variabil¬ 
ity. The p ower spectra of the LHS and hard-intermediate 
state fe.g.. lMendez fc van der Kli^llOOTl i lBelloni et al.ll200^ 1 
are dominated by a strong band-limited noise compo¬ 
nent that can r each fractional RMS amplitudes of up 


to ~ 50% (see l(V 

(1 6ti 3, 

.1 I 1971 I: iMendez & van der KlisI 

I 1997 I: IBelloni et al. 

1997: 

van der KlisI 20061. and refer- 


ences therein). Sometimes, low frequency quasi-periodic 
oscillations (QPOs) are observed with frequencies in the 
range of ^ 10~^ — 20 Hz. The characteristic frequencies 
llBelloni et al.ll2002h of thes e QPOs and noise components 
are found to correl ate (e.g., IWiinands fc van der Kliji 19991 : 
IPsaltis et al.l[l9^ : iBelloni et al.ll2002ll . generally increasing 


towards softer states. The soft-intermediate state shows no 


strong band-limited noise component (e.g.. 

Oda et alJll97(ll: 

IHoman et al.l I 2 OOII: IHoman & Bellonil 12005 

), and transient 


QPOs appear with frequ encies that are rather stable (e.g, 
ICasella et al.lf20o3 . 1200,^ . In the HSS, only weak power-law 


noise is observed in the power spectrum and sometimes a 
QPO with frequency in t he range 10-20 Hz (see, e.g., the 
HSS of XTE J1550-564; iHoman et ahlliooH l. Overall, the 
most common QPOs in BHCs have been classified based on 



Figure 1. Hardness-Intensity diagram of the four best sampled 
outbursts of the black hole GX 339—4. Hard colour is defined as 
the count rate in the 16.0—20.0 keV band divided by the count 
rate in the 2.0-6.0 keV band. Intensity is the count rate in the 
2.0-20 keV band. Both the hard colour and the intensity are nor¬ 
malized by the Crab (see Section[2]for more details). The four dif¬ 
ferent outbursts are plotted using different colours (see legend). 
Note that for the 2004-2005 outburst (red points) the hard-to- 
soft transition happened at a much lower (factor of ~ 5) intensity 
than for the other three outbursts. For the other three outbursts 
the transition occurred at approximately the same intensity. The 
rectangular shaded area shows the part of the outbursts that we 
used in the rest of the paper (see Figure [21l. The few gray points 
in this plot indicate observations in which the source was in the 
high-soft state, which we did not analyse in this paper (see Sec¬ 
tion [2]l. 


their characteristics as Type-A, B and C, and are u sually 
associated to different states (e.g.. ICasella et ^l2005h . 

One important aspect of the variability of any astro- 
physical signal is the possible phase lag between measure¬ 
ment s in two different ene rgy bands (e.g.. |Priedhorskv et al.l 
I 1979 I : iNowak et al.l Il999l ). The phase lag is a Fourier- 
frequency-dependent measure of the phase delay between 
two concurrent and correlated signals, in this case light 
curves of the same source, in two different energy bands. 
The phase-lag spectrum of the broad-band noise com¬ 
ponent i n the power density spectrum of black hole 
systems llMivamoto et akl Il98d iNowak et al.l Il999l) were 
first interpreted as due to soft disc photons that are 
Compton up-scattered in a corona of hot electrons that 
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Figure 2. Panel A shows a portion (see shaded area in Figure 0 of the full hardness-intensity diagram (HID) of the four best sampled 
outbursts of the black hole GX 339—4 with RXTE. Grey points indicate observations which were not studied in this paper (see Section!^. 
The black, red, green and violet points indicate the part of the four outburst (see text for details, and the legend for the identification of 
each outburst) for which we studied the lags. Hard colour is defined as the count rate in the 16.0-20.0 keV band divided by the count 
rate in the 2.0—6.0 keV band. Intensity is the count rate in the 2.0-20 keV band. Both the hard colour and the intensity are normalized 
by the Crab (see Section[2]and Figure^for more details) The hard-to-soft transition during the 2004-2005 (red) outburst took place at 
an intensity that is a factor ^^5 lower than in the other three outbursts. Panel B and C show the phase lags-intensity diagram (LID) for 
the HI (~ 5.7 — 15 keV) and H2 bands (~ 16 — 35 keV), respectively. The lags correspond to the average phase lag divided by 27r in the 
0.008 — 5 Hz frequency range using as reference the soft band (~ 2 — 5.7 keV). Horizontal arrows in Panel B mark the moment where the 
lags stop increasing with intensity (see Section [3] for more details). The vertical dashed-lines show that, as GX 339-4 becomes brighter 
in each outburst, the spectrum becomes softer (panel A) and the lags increase (panel B & C). The horizontal dotted-lines are shown to 
help the eye link the evolution of the lags as a function of the position of the source in the HID. In all panels, each point corresponds to 
the average value calculated from a full RXTE observation. 


surrounds the s ystem (see, e . g., Thorne fc Pric3 Il975l : 


IKazanas et al.l 119971. H oweve r, see, e.R., 
Now^i_rtalJ_(|l999|}j_lMa£cai;one_rtaL|_jl2000j); Pm^anerl 


Hua et all 199 j: 


(l200ll) : I Arevalo fc Uttlevl ~ 2006l ~) : lUttlev et alJ ( 201lK for 


problems with this interpretation). Phase lags from QPOs 
and the broad-band noise have been relatively well stud¬ 
ied in the LHS of a few sources, but using only a handful 


of observations (e.g 

van der Klis et al.l 19871: Nowak et al.l 

1999l:IWiinands et al. 

1 19991: iFord et al.ll 19991: Gui et al.ll200nl: 

Reig et al.l 20001: iMendez et al.l 2013tl. If the geometrv of the 


system changes significantly during an outburst, one would 
naturally expect that the phase lags between soft and hard 
photons will change as a function of the position of the 
source in the HID. To our knowledge, only the lags of the 
persistent (and somewhat special) BH system Cygnus X- 
1 have been studied in detail, taking advantage of the full 
Rossi X-ray Timing Explorer (RXTE) archive. As expected, 
the lags change as a f unction of the positio n of the source 
in the HID (see, e.g., iGrinberg et al.ll2014l . and references 
therein). 


In this paper, we present the first multi-outburst study 
of the time lags of the low-frequency variability for a tran¬ 
sient BHC LMXB. For this exploratory work, we chose 
the BHC GX 339-4, as evidence for lag evolution has 


been reported for the 2002/2003 outburst (see Figure 5 in 
iBelloni et al.l[2005h . GX 339-4 is p robably one of the best 


19971; 

Homan et al 

20051: Coriat et al. 

1200!^: Motta et al. 

20 111; 

iBuxton et al. 

I 2 OI 2 I: iRahoui et al 

I 2 OI 2 I: ICorbel et al. 


20 ll and references therein), and its outbursts are con- 


sidered as representative of the black hole population (e.g., 
IBelloni et al.fl200^ : Ivan der KlisI l2006l ~l . 


2 OBSERVATIONS AND DATA ANALYSIS 

We used all the 1414 pointed RXTE observations of GX 339- 
4 available in the HEASARC archive. In particular for this 
study, we use d data taken with the Proportional C ounter 
Array (PGA; IZhang et aI.lll9Si^ : Ijahoda et al.ll2006l l over a 
time span of ~14 years. To calculate X-ray colours and inten¬ 
sity, we used the 16-s time-reso lution Standard 2 mode data 
and the procedure described in lAltamirano et al.l (l2008ll : for 
each of the five PGA (Proportional Counter Units; PCU’s) 
detectors we calculated a hard colour defined as the count 
rate in the 16.0-20.0 keV band divided by the count rate in 
the 2.0-6.0 keV band. We defined the intensity as the count 
rate in the 2.0-20 keV band. Count rates in these exact en- 































































































































4 Diego Altamirano and Mariano Mendez 



Figure 3. HI phase-lag (top; lags are in the ~ 5.7 — 15 keV band respect to the ~ 2 — 5.7 keV band. See Section [2] for more details) 
and power spectra (bottom; 2—60 keV) of four representative observations. The red, green, blue and black colours mark observations 
92035-01-01-02, 92035-01-02-04, 92035-01-02-06 and 92428-01-04-00, respectively. The location of these observations in the HID and LID 
are indicated in Figure!^ 
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Figure 4. Hardness-intensity diagram (right) and HI Phase-Intensity diagram (left) for the 2006-2007 outburst. Hard colour is defined as 
the count rate in the 16.0-20.0 keV band divided by the count rate in the 2.0—6.0 keV band. Intensity is the count rate in the 2.0-20 keV 
band. Both the hard colour and the intensity are normalized by the Crab (see Section[2]for more details). The HI phase-lags (~ 5.7 — 15 
keV) are calculated relative to the (~ 2 — 5.7 keV) reference band, and correspond to the average phase lag divided by 27r in the 0.008 — 5 
Hz frequency range. The red, green, blue and black colours correspond to observations 92035-01-01-02, 92035-01-02-04, 92035-01-02-06 
and 92428-01-04-00, respectively. Their Phase spectrum and Power spectrum are shown in Figure [3] 
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Figure 5. 0.008 — 5 Hz frequency average phase-lag divided by 27r 
of band HI (~ 5.7—15 keV) relative to the ~ 2 —5.7keV soft band 
(coloured points; left y axis) and hard colour (gray points; right y 
axis) vs. the average 0.008 — 5 Hz fractional RMS amplitude (2-60 
keV) for the four outbursts discussed in this paper (see Figure [2] 
for the identification of the outbursts). Each point corresponds 
to the average values per observation, and for each measurement 
of a lag we plot the corresponding hard-colour point in gray. The 
hand-drawn dashed lines show the break in the hard-colour vs. 
RMS relation at RMS ~20%. The vertical and horizontal dashed 
lines are plotted to help the eye. 


ergy ranges were obtained by interpolating between PCU 
channels. We corrected the data for dead-time, subtracted 
the background contribution in each band using the stan¬ 
dard bright source background model for the PCA, version 
2 .lfl and removed instrumental drop- outs. In order to cor - 
rect for instrumental gain changes fsee l.Tahoda et aDl200fil l. 
we normalized all count rates by those of the Crab Nebula, 
obtained in the same energy range, same gain epoch and 
using the Crab observations closest in time to those of each 
GX 339-4 observation. For each observation we then cal¬ 
culated colours and intensity for each 16-s interval, and we 
averaged the colours and intensity over all PCUs. Finally, 
we averaged the 16-s colours per observation. 

For the timing analysis we used data from the PCA 
Event, Good Xenon and Single Bit modes. We calculated 
Leahy-normalized power density spectra (PDS) using all 
photons in the PCA bandpass for data segments of 128 sec¬ 
onds at the maximum time resolution available (except in 
the case of Good Xenon mode, where we binned the data to 
a resolution of 1/8192-s). No background or deadtime cor¬ 
rections were applied to the data prior to the calculation 
of the PDS. We subtracted a Poisson noise spectrum that 
includ es the effect of dead time based on the analytical func¬ 
tion of IZhang et al.l (Il995l l. and t hen converted the r esulting 
PDS to squared fractional RMS dvan der KliElll9^ 1. 

For each observation we also computed the complex 
Fourier transform separately for all photons in the absolute 
channel bands 0 — 13, 14 — 35 and 36 — 80; in the rest of this 
paper we refer to these bands as soft, HI and H2, respec¬ 
tively. These bands correspond to approximately the energy 
ranges 2 — 5.7 keV, 5.7 — 15 keV and 15 — 34 keV, respec¬ 
tively, and were selected in order to maximize the number 
of observations where we could use the exact same chan¬ 
nel bands. While we used always the same channel-ranges, 
the corresponding energy-ranges v aried slightly in tim e due 
to the instrumental gain changes djahoda et al.ll2006fl . We 
estimated that those variations affect the lags by less than 
~ 10%. We calculated average cross spectra using the soft 
band as the referen c e ban d, an d following th e descr iption in 
I Vaughan fc NowakI dl997h and iNowak et all dl999l L we cal¬ 
culated the phase lags as a function of Fourier frequency for 
each of the cross spectra. We use the nomenclature HI and 
H2 as well to refer to the lags between photons in the HI 
or H2 bands, respectively, relative to the photons in the soft 
band. Given our convention, positive lags means that hard 
X-ray photons lag the soft ones. 

A visual inspection of our data indicates that the shape 
of the phase-lag spectra is complex and dependent on the 
observation. One could divide each lag spectrum into dif¬ 
ferent components, as it is commonly done for the power 
spectra of these same observatio ns (generally in the form of 
mult i-Lorentzian model, see, e.g., lNowakl2000l : [Belloni et m] 
I2OO2II . However, the quality of the lag-spectra does not al¬ 
ways allow to differentiate between components; therefore, 
as a first approach, in this paper we study the evolution of 
the frequency-averaged 0.008-5 Hz phase-la£|. This type of 

^ PCA Digest at http://heasarc.gsfc.nasa.gov/ for details of the 
model 

^ We calculated the frequency-averaged 0.008-5 Hz phase-lag, 
by first averaging the raw cross-correlation vectors, then aver¬ 
aging over frequency, and finally calculating the lag as PL = 
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analysis is analogous to those works where the relation be¬ 
tween average fractional rms amplitude of the power spectra 
of BHs (and also NS) is studied as a function of the spectral 
state of a sou rce (e .g. |Belloni_gt__alJ|2()0^; Fender^^^ij|l200^ 
Motta et al. ID; iMufioz-Darias et alJ l^llal : IStiele et all 


201 ih . A possible multi-component description of the phase- 


lag spectrum will be discussed in a separate paper. 

Dead-time driven cross talk between different energy 
channels induces a ±7r phase lag in the Poisson -noise domi¬ 
nated part of the PDS (Ivan der Klis et aDllOSTl l. To correct 
for this, we subtracted the cross spectrum in a frequency 
range well above that of the broad-band noise (200 — 1000 
Hz), in which the X-ray variability is dominated by Pois¬ 
son noise. We inspected and confirmed that, as expected, 
the cross spectrum in this frequency range has no signif¬ 
icant imaginary part; hence this procedure has no effect 
on the sign of the phase lags (see Ivan der Klis et ^ Il987l . 
for more details). Unless stated explicitly, quoted errors use 
= 1.0. 


Throughout the paper we give phase lags in units of 
rad/27r (cycles) and we do not report on time-lags (defined 
as phase-lag divided by the frequency u), as the frequency 
range we use (0.008-5 Hz) spans nearly a factor ~600 so the 
average time-lag will strongly depend on the binning used 
(e.g. linear or logarithmic), and can be strongly biased by 
the lags at frequencies lower than 1 Hz. The average phase- 
lags measurements we report here can be used to calculate 
time-lags at a given frequency to a first-order approximation. 


3 RESULTS 

We were able to measure the phase-lags in the 3 bands de¬ 
fined in Section[2]in 1362 out of the 1414 observations avail¬ 
able in the archive. Approximately 1020 out of these 1362 
observations sample four complete outbursts of GX 339- 

4 that showed the canonical full loop in the HID; the 
HIDs of the four outbursts studied here are plotted in Fig¬ 
ure [p As described in Section [J (and in previous works, 
e.g.. lBelloni et ^l2005l : iMotta al.ll201ll . and references 
therein) and shown in Figure [T] the outbursts of GX 339-4 
evolve in an anti-clockwise “q-shaped” form. 

In Panel A of Figure[2]we show a portion of the HID of 
these four outbursts. For each outburst the gray points show 
the full outburst, while we indicate the rising part of each 
outburst (hard colour > 0.2) with the black, red, green and 
violet points for the 2002-2003, 2004-2005, 2006-2007, an 
2010-2011 outbursts, respectively. The transitions between 
the hard and the soft state are well sampled with RXTE 
observations, except for the 2002-2003 outburst (black), in 
which the hard-to-soft transition was sampled by only 5-6 
observations. In panels B and C of Figure [2] we plot the 
source intensity versus the 0.008-5 Hz phase-lags for the HI 
and H2 bands, respectively; in both cases we use the soft 
band as reference. We will refer to these plots as the (Phase- 
)Lag-Intensity Diagrams, LID. 

For the 3 brightest outbursts. Panel B in Figure [2] shows 
that the evolution of the HI lags can be divided into 3 inter¬ 
vals: initially the lags show a slow, but signihcant, increase 


arctan (3/3f{)/27r, where 3 and 5R are the imaginary and real part 
of the resulting frequency-averaged vector, respectively. 


with intensity until the source “turns the corner” in the HID; 
at that point the lags start to increase at a much higher rate 
with intensity. Finally, the lags appear to reach a maximum 
at -^0.03 cycles , decrease to 0.020 — 0.025 cycles, and remain 
approximately constant with intensity. After this point we 
can no longer measure the lags in single observations, or the 
source starts undergoing fast state transitions. For the sake 
of clarity, this more complex behaviour will be presented in 
a followup paper. 

For the 2002-2003 (black), 2006-2007 (green) and 2010- 
2011 (violet) outbursts, the HI lags reach the maximum 
value on MJD 52400-52402 (ObsID: 70109-01-06-00/10-00), 
MJD 54140.2 (ObsID: 92035-01-03-00 and MJD 55298.7 
(ObsID: 95409-01-14-03), respectively (marked with hori¬ 
zontal arrows in Figure[2l Panel B). These dates correspond 
to the moment when the Type-C QPO reaches frequencies 
in the 2 — 3 Hz range, close to, but still below, our 5-Hz 
limit for the calculation of the lags. To investigate whether 
this could be the reason for the behaviour of the lags at this 
point in the outburst, we also calculated the average phase- 
lag in the 0.008 — 10 Hz range (not shown); in this case we 
found that the overall pattern in the LIDs is the same as 
those shown in Panel B in Figure [2] with a maximum of 
~0.03 cycles happening at the same position of the arrows 
in Panel B. In none of the observations of the 3 brightest 
outbursts we detected a Type-B QPO; these QPOs occur 
when the source softens fur ther in the outburst (see, e.g., 
table 3 in iMotta et aklEoilh . 

In the case of the 2004-2005 outburst (red), the lags 
increase with intensity during the whole rising part of the 
outburst. The “turn the corner” point occurs at an intensity 
of ~0.12 Grab, and the softening in the HID is less sharp 
than in the other 3 outbursts. This might be connected to 
the fact that this is also the only full outburst which made 
the hard-to-soft transition at a significantly lower luminosity 
than the others. In this outburst the lags in the HI band are 
the longest (up-to ~0.045 cycles). The last two observations 
presented here (ObsID:90110-02-01-03/02-00), the one with 
the longest lags, already contain Type-B QPOs in the power 
spectrum. 

The evolution of the H2 lags during the individual out¬ 
bursts in Panel C is similar to what we observe in Panel 
B, with the main differences being that: (i) the lags dur¬ 
ing the 2004-2005 outburst (red) show a maximum of ~0.17 
cycles (at MJD 53232.35, ObsID:90110-02-01-00/02), and a 
sudden drop after that; the two observations after the max¬ 
imum in this Panel correspond to those where we detect 
Type-B QPOs (see above), and (ii) the lags measured in the 
other 3 outbursts show a fast increase at the very end of 
our sampling, reaching a maximum in the ~0.14-0.17 cycles 
range. 

In all outbursts, the point where the lags start increas¬ 
ing at a higher rate (see panels B & C in Figure [2]) coincides 
with the “turn the corner” point in the HID (see panel A), 
where the source starts the rapid transition to the inter¬ 
mediate state (in all panels of Figure [2] we plot horizontal 
dotted-lines to guide the eye), strongly suggesting that both 
phenomena are connected, and related to the same changes 
the system is undergoing during this softening. This “turn 
the corner” point in the HID does not appear to be related 
to the appearance of the Type-C QPO, as at least in the 
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2006-2007 and 2010-2011, clear Ty pe-C QPOs a re already 
detected before this transition fe.g.. iMotta et aDl201ll '). 

In Figure [5] we plot the HI lags (left y axis) vs. the 
fractional RMS amplitude computed in the 0.008 — 5 Hz 
range, using the full 2-60 keV RXTE energy band (coloured 
points). To compare with the correlatio n between fractiona l 
RMS amplitude and hard colour le.g.. iBelloni et ^l2005ll . 
we also plot in this Figure the relation between hard colour 
(right y axis) vs. fractional RMS amplitude (gray points). 
For the 3 brightest outbursts, the point where the HI lags 
reach their maxima (panel B of Figure[2]) corresponds to the 
point where the average RMS amplitude is ~ 20%. Again, 
this is not the case for the 2004-2005 (red) outburst, where 
the HI lags continue increasing up to the end of our sam¬ 
pling. From Figure[S]it is also apparent that the hard-colour 
vs. RMS relation has a break when the RMS amplitude is 
>^20% (see hand-drawn dashed lines in each panel), even 
for the 2004-2005 outburst. (This break was already appar- 
en t, but not specifically discussed, in some of the Figures 


Belloni et alJl2005l : Il4nder et al.(l2009l : iMotta et al.ll201l] : 


IStiele et al.ll201lfr 


We not e that GX 339-4 h as undergone additional out¬ 
bursts fe.g. IMotta et al']l201ll ). however we do not report 
them here in detail as they were either outbursts where 
GX 339-4 always remained in the hard state, or the outburst 
was not well sampled by RXTE observations. Nevertheless, 
the lags during the hard state of those other outbursts are 
consistent with the ones shown in Figure[21 In Figure[21 when 
hard colour < 0.2, the RMS amplitude in the PDS decreases 
further (not shown in Figure [5] see e.g.. IMotta et al1l201ll ~l 
and, as mentioned above, the lags in single observations are 
not always well constrained. During the soft-to-hard tran¬ 
sition at the end of these outbursts, t he RMS amplitude 
in the PDS incr eases once more (e.g., IBelloni et al.l I^OSl : 
IStiele et al.ll201lll . and we can again measure the lags in in¬ 
dividual observations; in this phase of the outburst the lags 
are consistent with those in the LHS during the rise of the 
outburst. We do not discuss those lags further in this paper. 


4 DISCUSSION 

We present the first systematic study of the evolution of 
the X-ray lags of the broad-band variability component 
(0.008 — 5 Hz) in the black-hole candidate GX 339—4 as 
a function of the position of the source in the hardness- 
intensity diagram (HID). The hard photons always lag the 
soft ones, with the phase lags ranging from ~ 0.005 to ~ 0.17 
cycles depending on the observation and the energy band 
used to measure the lags. In all four outbursts studied here, 
as the source brightens at the beginning of the outburst, the 
lags initially increase slowly as the source is in the low-hard 
state (LHS), and then they start to increase much faster with 
intensity as the source initiates the transition to the hard- 
intermediate state (HIMS). After reaching a local maximum, 
the lags decrease and remain more or less constant (Fig¬ 
ure m, or decrease and then increase again (Figure [5p), 
before the source moves into the high-soft state (HSS). While 
there is no unusual feature in the HID at the time that the 
lags reach the maximum value, the maximum of the lags 
coincide with a signihcant break in the relation between the 
hard-colour and the RMS amplitude of the broad-band com¬ 


ponent in the power spectrum of the source. For the 2002- 
2003 outburst, this point also coincides with the time at 
which liioman et al.l (I 2 OO 5 I I observed a dramatic drop of the 
optical/NIR flux which they interpreted as a sudden change 
in the properties of the compact jet in this source. 


4.1 Phase-Lags and the Optical/NIR behaviour of 
GX 339-4 

Several works have studied the relation between the X- 
ray emission and th at at other wavel e ngths in Blac k Hole 
binaries (see , e.g. , Ijain et H] I 2001 I : iRussell et al.l I 2 OIII : 
Buxton_et_^lJ[20li and references therein). In particular, 
Homan et al.l (I 2 OO 5 I I studied the X-ray and Optical/Near- 
Infrared evolution of the 2002-2003 outburst of GX 339-4 
(black points in our Figures[2]andO. They found that in the 
LHS (MJD 52330-52398), both the optical/NIR and X-ray 
fluxes increase and remain well correlated with each other. 
As GX 339-4 evolved into the intermediate state, starting on 
MJD 52400, the optical/NIR fluxes decreased rapidly until 
GX 339-4 reached the spectrally soft state (about 10 days 
later), where the opti cal/NIR fluxes rema ined low and ap¬ 
proximately constant. lHoman et al. I ll2005l) noticed that dur¬ 
ing the intermediate state, and particularly on MJD 52402, 
the disc flux increased dramatically (from < 1% to about 
30% of the total flux; note however that the total flux in the 
RXTE band actually decreased slightly at this point) with¬ 
out substantially changing nor interrupting the evolution of 
the frequency of the Type-G QPO or the power-law index 
in the X-ray spectra. Unfortunately there is a 3-day gap, 
between MJD 52402 and 52405, in the RXTE coverage of 
the 2002-2003 outburst (black points in Figures [2] and [5}. If 
the lags in this outburst evolved in the same way as in the 
2006-2007 (green) outburst, where the HI lags drop from a 
local maximum of ~0.03 cycles to ~0.02 cycles within one 
day, we would have probably observed a maximum in the 
HI lags in panel B (and possibly also in the H2 lags, panel 
G) within a day or two of MJD 52402. If that was the case, 
this would indicate that in GX 339-4 there is a relation be¬ 
tween the lags and the optical/NIR emission, whereas the 
latter is not relate d to the power-law index or frequency of 
the Type-G QPO llHoman et ahllioosh . The RMS amplitude 
of the variability would also be related to the optical/NIR 
emission, as there is a clear break at 20% RMS amplitude in 
the hard colour vs. RMS relation, coincident with the point 
where we observe the maximum in the HI lags (Figure [S] 
see below). 

The accelerated steepening of the power-law compo¬ 
nent in the X-ray spectrum around MJD 52400, together 
with GX 3 39-4 showing a dram atic decrease in Optical/NIR 
fluxes, led iHoman et al.l ll2005ll to suggest that this is the 
moment where the characteristics of t he compact jet (as 
usually observed in the hard state, e.g.. lGorbel et aDl2013fl 
start to change. If the corona of hot Gomptonizing elec¬ 
trons and the compact radio jet are related, or if the base of 
the jet is the p l ace where Gomptoniz a tion takes place (e.g ., 
iGiannios et al.l l2004l : iGiannio a I 2 OO 5 I : iMarkoff et al.l 12003 1 , 
one would expect that the properties of the jet and the lags 
are correlated. Indeed, given that in BHCs in the LHS the 
radio and X-ray fl ux are correlated (also for GX 339-4, e.g., 
iGorbel et al.ll20f^ . and references therein), and that the lags 
correlate with X-ray intensity (Panel B & G in Figure [2]), 
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Figure 6. HI phase-lag (~ 5.7— 15 keV band respect to the ~ 2 — 5.7 keV band) vs. hard colour (normalized by the Crab) for Cyg X-1 
(gray) and GX 339-4. Hard colour is defined as the count rate in the 16.0-20.0 keV band divided by the count rate in the 2.0—6.0 keV 
band. 


indicates that, at least in the LHS, the lags probably corre¬ 
late with the radio flux of the compact jet. Several spectral 
and timing properties of accreting systems correlate with 
one another; here we have shown that the phase lags also 
correlate with those properties. To progress further, and to 
gain insight on what are the basic source properties that 
drive these correlations, one would need to test, for instance, 
whether the lags in GX 339-4 (and other sources) also cor¬ 
relate with the radio flux in the intermediate state, where 
the lags do not (always) correlate with the source intensity 
(see, e.g., Figure [2]). Such a study could prove whether the 
lags as we estimated in this work are a good indicator of 
some of the radio jet characteristics, similar to the evidence 
provided by the correlations between radio luminosity and 
X-r ay timing frequencie s found in both NSs and BHs (see, 
e.g., Migliari et akllioosh . 


4.2 Lags and the fractional RMS amplitude 

In Figure [5] we showed that there is a break in the hard- 
colour vs. RMS correlation that occurs consistently at RMS 
~ 20% in all four outbursts. This break can also be seen as 
a “notc h” in the absolute RMS-inte nsity diagram (e.g.. Fig¬ 
ure 1 in lMufioz-Darias et al.l[2011al . at absolute RMS ~ 220 
cts/s and PCU2 count rate ~ 1000 cts /s, corresponding to 
a fractional RMS amplitude of ~ 20%). Munoz-Darias et al.l 


ll2011ah show that, at least in the 2006-2007 outburst, 
this notch in the absolute RMS-intensity diagram coin¬ 
cides with the observations where X-r ay spectral fits re ¬ 
quire a thermal blackbody component (iMotta et al.ll201ll ): 
iMunoz-Darias et al.l ll2011ah interpreted this as the appear¬ 
ance (in the 2-15 keV band) of an optically thick accre¬ 
tion disc with a very low variability level. For the 2002-2003 
(black) outburst, the break in the hard-colour vs. RMS rela¬ 
tion also coincides with the time when the disc flux showed 
a dramatic increase (from < 1% to about 30% of the total 
flux), and the optical and NIR bands starte d to decrease 
(see discussion above, and iHoman et al.l[2005l ). As it can be 
seen in Figure [5] except in the case of the weak 2004-2005 
(red) outburst (see Section[3]and Figure[5]for details), in the 
other three outbursts the break in the hard-colour vs. RMS 
relation coincides with the moment when the HI lags reach 
their maximum, again suggesting a strong relation between 
lags, RMS amplitude, and spectral evolution. 

As it is apparent from Figure 5 in iFender et al.l ll2009l) . 
at least two other BHC systems (XTE J1550-564 and XTE 
J1859-I-226) also show a break in the hard-colour vs. RMS 
relation that, coincidentally, happens at around 20% RMS 
as well. Although this coincidence has to be explored in 
more detail (as th e frequency and energy bands used by 
iFender et al.|[200^ . differ from those used in this paper), it 
suggests the interesting possibility that we can identify the 
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MJD - 50946 (days) 


Figure 7. Cut-off energy Ecu t of the energy spectra (top) as 
reported bv iMotta et alj i2009^ and HI phase-lags (bottom; lags 
are in the ~ 5.7 — 15 keV band respect to the ~ 2 — 5.7 keV 
band. See Section [2] for more details) during the first 45 days of 
the 2006-2007 outburst of GX 339-4. With continuous lines we 
show a broken power law fit done to both curves simultaneously, 
with the condition that break is the same for both fits and that 
the power-law index before the break is zero. For this model, the 
break is at MJD —5094,6 = 20.1 ±1 day (see Section liTil for more 
details). 


point where the spectral fits require a thermal component 
just by the use of the fractional RMS amplitude or the lags, 
assuming that other sources show the same evolution of the 
lags during the outburst. 


4.3 GX 339—4 vs. Cygnus X—1 

iGrinberg et al.l (l2014l ) reported on the long term evolution 
of the energy-resolved X-ray variability of the persistent 
BHC LMXB Cygnus X-1 using data taken between 1999 
and 2011. These authors find positive time lags that evolve 
as a function of the spectral state of the source: the av¬ 
erage time lags (as calculated in narrow frequency ranges) 
increase considerably when the source moves from the hard 
to the intermediate state (with a maximum of ~20 msec 
for the 0.1-30 Hz lags of the 9.4-15 keV photons relative 
to the 2.1-4.5 keV p hotons; bands 4 and 1, respectively in 
ICrinberg et al.ll2014h . and then the lags suddenly drop to 
< 5 msec when Cygnus X-1 en ters the softest states . To 
compare our results with those of ICrinberg et aP ll2014h . we 
calculated the average 0.008-5 Hz phase-lag for Cyg X-1 us¬ 
ing the same energy bands and the exact procedure used 
for CX 339-4. In Figure [6] we show the phase Lags of Cyg 
X-1 (gray point^ and CX 339-4 (colour points; colours are 
the same as in Figure [SJ. This Figure shows that Cyg X-1 
and CX 339-4 span the same range in hard colour (although 
note that CX 339-4 soft-state data are not included in this 
plot). The overall shape drawn by Cyg X-1 and CX 339- 
4 is similar: As the source softens, the lags increase until 


® Note tha t the datapoints form a pattern very similar to those 
showed by IGrinberg et al.l ||2014|1 , although we are plotting an 
averaged phase-lag i n the 0.008-5 Hz frequency range, while 
IGrinberg et al.l ll2014^ were plotting average time-lag in narrow 
frequency ranges. 


a maximum is reached, and the lags decrease. In Cyg X-1 
the turn over happens at softer colours than in CX 339-4, 
and the decrease in lags continue to much lower values than 
in CX 339-4. The differences in absolute value of phase-lag 
(and that of the hard color at which the lags change sharply) 
could be related to differences in the system’s characteris¬ 
tics (e.g. black-hole mass) or to the fact that Cyg X-1 does 
not r each the typical soft state seen in other BH systems 
(e.g., ICrinberg et al.ll2014l . and references therein). In any 
case, the fact that Cyg X-1 shows the turn-over as it tran- 
sits to a softer state while it does not show any QPOs (see 
ICrinberg et al.ir2014l . for a discussion) supports our conclu¬ 
sion that the turn-over in phase lags observed in CX 339-4 
is not directly related to the appearance or disappearance 
of QPOs (see Section O. It remains to be seen whether the 
lags also drop in CX 339-4 when it reaches the soft state. 


4.4 On the origin of the lags 


In CX 339-4 the photons in both hard bands, HI and 
H2, lag the photons in the soft band. This is consistent 
with previous r esults of the lags in the LHS of this and 


other BHC (e . g -. vandn£^llset5][l987j^ Nowak_et_al 


w 

Cui et al.ll2000l: Munoz-Daxias et aLll2011bl : ICrinberg et ^ 


Belloni et al. 2005: Hua et al.l 19971: jKazanas e^ 


1999 


1997 


2014h and NS iFord et al.lll999l) LMXBs . 

Hard lags have been originally explained as due to in¬ 
verse Compton scattering of soft disc photons in a uni¬ 
form cloud of hot electrons (the so-called corona) close 


to the compact object 


Pwne| 1980|; iMivamoto et al.l 


ll988l : l^ord et al.lll999l : iNowak et al.lll99^ ). With each scat¬ 
tering the energy of a photon {E) i ncreases approximately 
by /S.E/E = kTe/{meC^) (e.g., IMivamoto et al.l 1198^ 
iNowak et aLlll99^ , where kT^ is the temperature of the elec¬ 
trons in the corona. This implies that the longer the path 
length of the photon has to travel before leaving the corona, 
the longer the time lag, and the higher the energy of the pho¬ 
ton. In this scenario the minimum time lag between hard and 
soft X-rays would be of the ord er of the light crossing time of 
the Comptonising re g ion (e.g. .IM ivamot o fc Kitamotd[l98^ : 
iKazanas et al.l 1 19971 : INowak et al.l 19991) . while the maxi¬ 
mum time lag must be less than the size of the emitting 
region divided by the slowest propagation speed (sound 
or th ermal wave speed) in that region (e.g INowak et al.l 
I 1999 I) . Very roughly, the average phase-lags in CX 339-4 
would imply a light-crossing time of ^21000 km, or ~700 
Schwarzschild radii {Rs) for a 10 Mq black hole (or much 
large r radii if one only takes the longest la gs observed, see, 
e.g., IKazanas et aL|[l9^ : iHua et al.^ll997^ . 


_ T he main obstacl e to this scenario ( e.g. Npwak__eLaL 

199 ^; iLin et aP I2OOOI : iMaccarone et all I2OOOI : iPoutaneii 


20011) is that in a constant-density corona the energy re¬ 


quired to maintain a population of highly energetic elec¬ 
trons is only available very close to the black hole, whereas 
the magnitude of the lags imply either very energetic elec¬ 
trons far away from the black hole, or a corona with an 
optical depth much larger than deduced from spectral fits 
llWardzihski et ahl I 2 OO 2 I : iMivakawa et al.l I 2 OO 8 I ) . The prob¬ 
lem becomes worse if, to explain the frequency dependence of 
the time lags in these systems, the density of the corona goes 
as r~^ (i.e. inversely proportional to the distance to the com - 
pact object, see, e.g.. IKazanas et al.lll997| : lHua et al. IH). 
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0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 

2-20 keV Intensity (Crab) 

Figure 8. The ratio between the lags in the H2 (~ 16 — 35 keV) and HI (~ 5.7 — 15 keV) bands as a function of Intensity (in units of 
the Crab Nebula). Lags correspond to the average phase lag divided by 27r in the 0.008 — 5 Hz frequency range using as reference the 
soft band (~ 2 — 5.7 keV). We only show the ratios for intensities higher than 0.1 Crab, as below that the lags in single observations are 
not well constrained (however, the average ratio between lags is still ~2). 


In addition. iMaccarone et ciD (l2000l ) showed that the energy 
dependence of the width of the autocorrelation function of 
the X-ray light curves produced by Compton scattering in 
an isotropic corona is inconsistent with the observations. 


_ These difficulties lent support to the idea of iLvubarskiil 

lll997h that the variability (and the lags) in the light 
curve of these systems could be due to small variations 
of the viscosity parameter in the accretion disc. These 
variations render fluctuations of the mass accretion rate 
that propagate through the disc, and eventually dissi¬ 
pate in the corona. Several authors explored the effect 
of t his idea on the lag spectrum of accreting sources 
Misr^ 2000l : iKotov et al.l 200lj^ Arevalo fc Uttlevll2006l : 
lUtUev et al.l 201ll) . Arevalo fc Uttle^ ( 20061 ) carried out a 
detailed study of the dependence of the power density and 
lag s pectrum (and coherence function; IVauehan fc NowakI 
[ 1993 ) upon, among others, the emissivity profile of the disc, 
and the geometry and the vi scosity parameter o f the accre- 
tion flow. Using this model, lArevalo fc Uttle"^ 1 I 2 OO 6 I) were 
able to reproduce the observed lag spectra in an RXTE ob¬ 
servation of Cyg X-1 and an XMM-Newton observation of 
the narrow line Seyfert 1 galaxy NGC 4051. However, while 
the same model can reproduce the 2 — 3 vs. 0.5 —0.9 keV lags 
observed with XMM-Newton in GX 339-4 in the low-hard 


state during the 2004 outburst, it fails to explain the 6 — 9 
vs. 2 — 3 keV in the same observation llUttlev et all gOffijb 
iGian nios et ^ (|200 4), bas ed on the work of Reig et ^ 
1 I 2 OO 3 I) and Kvlafis et al. I 2008h . proposed that the hard lags 
could be due to Comptonisation of soft disc photons in the 
jet that is present in the low-hard state of galactic BHGs. 
The proposed mechanism is the same as in the case of Comp¬ 
ton scattering in a corona (see above), but in this case the 
anisotropy of the scattering process along the jet resolves 
the energy problem compared to the case of the lags pro¬ 
duced in a (spherical) coronc0. Furthermore, in this model 
the width of the autocorrelation functio n decreases with 
energy, in agreement with observations (IMaccarone et al] 


lUttlev et al.l (l201lh argued that, contrary to what they ob¬ 
served in GX 339—4, this model predicts that the 2 — 3 keV vs. 
0.5 — 0.9 keV lags should be smaller than the 6 — 9 keV vs. 2 — 3 
keV lags due to dilution by the direct seed photons from the disc 
at soft energies. However. [Kara et al.l 1 I 2 OI 3 II recently showed that 
contamination by either Poisson noise or an incoherent compo- 
nent does no t affec t the lags of a coherent signal. The findings of 
lUttlev et al.l ll201lh would still be consistent with the jet model 
of the lags if, for instance, the fluctuations in the disc were Pois- 
sonian, while the observed variability was due to the response 
(Green) function of the Comptonising component (the jet in this 
case). 
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I 2000 I) . It is interesting that the radio emission from the 
jet is quenched when the source m oves from the low- hard 
to the soft-intermediate state ('e.g.. Ihhnder et ^l2009l . and 
references therein), while coincidentally, in GX 339-4 we 
find that the lags drop as the source moves into the hard- 
intermediate and the soft-intermediate state. 

The sudden change we observed in GX 339-4 is sim¬ 
ilar to that seen in Gyg X-1, for which the magnitude of 
the lags drop abruptly as the source moves out of the hard 
state (see Figure [6]). In the model where the lags are pro¬ 
duced by Comptonisation in the jet, the increase of the lags 
as the intensity of the source increases in the LHS of GX 
339-4 (see Figure [2]) would be due to an increase of the 
he ight and/or the radiu s of the base of the jet (see figure 5 
in iGiannios et al]l2n04h . On the other hand, the drop of the 
lags could indicate that the optical depth of the Comptonis- 
ing region in the jet drops, consistent with the interpretation 
that the spectrum of the radio jet changes from optically 
thick to optically thin sy nchrotron emission at this point 
in an outburst (see, e.g., iFender et al.l 1200^ : ICorbel et al.l 

I 2 OI 3 II . 

Interestingly, the maximum of the lags in GX 339-4 
coincides, within ±1 day, with the time at which there is 
a sudden incre ase of the cut-off e nergy Ecut of the cut-off 
power law that iMotta et al] (l2009ll fitted to the energy spec¬ 
trum of GX 339-4 during the 2006-200 7 outburst (se e their 
figure 6). Furthermore, from figure 6 in iMotta et al.l ll2009l l 
it is also apparent that the Ecut of the hard spectral com¬ 
ponent in GX 339—4 remained more or less constant for the 
first ~ 20 days of the 2006-20076 outburst, after which, Ecut 
started to decrease. From our Figure U it appears that the 
behaviour of the phase lags also changes around that time; 
at the beginning of the outburst the phase lags appear to 
be constant, and when the intensity reaches ~ 0.2 Grab the 
lags start to increase . 

To further investigate this, in F igure [3 we plot the tem¬ 
poral evolution of Ecut (data from iMotta et al.l l2009l l and 
the phase lags for the first 45 days of the 2006-2007 out¬ 
burst. This Figure shows that both Ecut and the phase 
lag remain more or less constant until approximately MJD 
50952, date when Ecut starts to decrease while the phase lags 
start to increase. Under the assumption that both quanti¬ 
ties are correlated, we fitted a broken power law to both 
curves simultaneously, with the condition that the break 
was the same for both fits; we further assumed that the 
power-law index before the break was zero, i.e. both Ecut 
and the phase lags are constant before the break. The con¬ 
tinuous line in both panels of Figure [7] represents the best¬ 
fitting model ix^/d-of = 4.1 for 30 dof) with a break at 
MJD —50946 = 20.1±1 day. The fit improves /dof = 2.6 
for 28 dof) if we let the power-law indices before the break 
free to vary during the fits; in this case the break happens 
at MJD — 50946 = 17.4 ± 0.7 day. None of the fits is sta¬ 
tistically acceptable and there is a priori no reason to as¬ 
sume that a broken power-law is the correct model to de¬ 
scribe the data. Even if our proposed model was correct, 
the lack of data at the break precludes us from conclud¬ 
ing whether the break occurs simultaneously in both curves. 
While this needs to be investigated in more detail (e.g., with 
data from other outbursts), Figure [7] strongly suggests that, 
already at this early phase of the outburst, both Ecut and 
the phase lags are correlated. This result, together with the 


fact that the maximum lag we measured in GX 339-4 coin¬ 
cides (within ±1 day) with the time at which there is another 
(sudden) increase of Ecut, suggests that the average magni¬ 
tude of the lags are related to the properties of corona/jet 
fe.g.. iG iannios et al.| 2004 ), rather than to the disc (e.g., 
iLvubarskiill 19971 : lUttlev et al.ll201l| j. 

Intriguingly, there is no apparent discontinuity in the 
time evolution of Ecut at the point at which the lags start 
to increase very rapidly, e.g. the green point during the 2006- 
2007 outburst in the left panel of Figure^] This point, how¬ 
ever, coincides with the turn-the-corner point in the HID 
(right panel of Figure 21 see Figure [7] for the other out¬ 
bursts), which defines the transition from the hard to the 
intermediate state. To investigate this in more detail, in 
Figure [S] we plot the ratio between the 0.008-5 Hz aver¬ 
aged lags in GX 339-4 measured in the H2 and HI bands 
as a function of the intensity of the source. At intensities 
lower than about 0.35 Grab (for the weak 2004-2005 out¬ 
burst this value is ~ 0.12 Crab) the ratio of the lags re¬ 
mains approximately constant at a value of ~ 2, and as 
the intensity rises further the ratio increases abruptly up 
to values of 6 — 7. The ratio in Figure [8] deviates from a 
constant value at the same time that the source reaches the 
“turn the corner” point in the HID (Figure [T]), implying 
that the transition to the intermediate state has an effect 
on the characteristic of the component in the accretion flow 
that prod uces the la g s, e.g ., the geometry of the jet in the 
model of iReig et ^ (1200^ and IGiannios et al.l ^2004!) , or 
the accretion disc in th e model of LYubarskii (Il997 1 and 
lArevalo fc Uttle"^ 1 I 2 OO 6 II . IMotta et al.l ( ^09l l found that, 
at least in the 2006-2007 outbursts of GX 339-4, the ac¬ 
cretion disc does not contribute significantly to the emis- 


sion in the RXTE/ PCA band in this phase of the outburst 
llMotta et al.ll2003 . used the same data that we used in our 
paper for the study of the lags). Together with our findings, 
this implies that, even if the lags are due to the propagation 
of fluctuations in a relatively cool disc that peaks outside the 
RXE/PCA band, Comptonisation would still be required to 
explain lags. 


5 SUMMARY 

We studied the evolution of the X-ray lags of the transient 
BHC-LMXB GX 339-4 as a function of the position of the 
source in the HID. In the low-hard state the lags correlate 
with X-ray intensity, and as the source starts its transition 
to the intermediate/soft states, the lags first increase faster, 
and then appear to reach a maximum, although the exact 
evolution depends on the outburst and the energy band used 
to calculate the lags. Changes in the lags appear to coincide 
with changes in Optical/NIR flux, fractional RMS ampli¬ 
tude, and probably the appearance of a thermal component 
in the X-ray spectra, strongly suggesting that lags can be 
very useful to understand the physical changes that GX 339- 
4 undergoes during an outburst. We find evidence for a con¬ 
nection between the evolution of the cut-off energy Ecut (of 
the energy spectra) and the phase lags, suggesting that the 
average magnitude of the lags are related to the properties 
of corona/jet, rather than to the disc. 

The behaviour of the lags in GX 339-4 i s sim ilar to 
that of the BH Cygnus X-1 dCrinberg et al.l 120141 '). sug- 
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gesting similar phenomena could be observable in other BH 
systems. Understanding how lags evolve in other BH out¬ 
bursts, and whether the correlation with intensity, maximum 
value at ~20% fractional RMS amplitude, etc., is common 
to all BHC outbursts, or whether that depends on the char¬ 
acteristics of the system (e.g., neutron star vs. black hole, 
BHC mass, spin, inclination of the orbit with respect to the 
observer, radio flux, etc), is a necessary input for current 
and future models that describe both X-ray variability and 
multi-wavelength observations. 
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